Analysis of signalling pathways using continuous time Markov chains

We describe a quantitative modelling and analysis approach for signal transduction networks. We illustrate the approach with an example, the RKIP inhibited ERK pathway [CSK+03]. Our models are high level descriptions of continuous time Markov chains: proteins are modelled by synchronous processes and reactions by transitions. Concentrations are modelled by discrete, abstract quantities. The main advantage of our approach is that using a (continuous time) stochastic logic and the PRISM model checker, we can perform quantitative analysis such as what is the probability that if a concentration reaches a certain level, it will remain at that level thereafter? or how does varying a given reaction rate affect that probability? We also perform standard simulations and compare our results with a traditional ordinary differential equation model. An interesting result is that for the example pathway, only a small number of discrete data values is required to render the simulations practically indistinguishable

Trend-based analysis of a population model of the AKAP scaffold protein

We formalise a continuous-time Markov chain with multi-dimensional discrete state space model of the AKAP scaffold protein as a crosstalk mediator between two biochemical signalling pathways. The analysis by temporal properties of the AKAP model requires reasoning about whether the counts of individuals of the same type (species) are increasing or decreasing. For this purpose we propose the concept of stochastic trends based on formulating the probabilities of transitions that increase (resp. decrease) the counts of individuals of the same type, and express these probabilities as formulae such that the state space of the model is not altered. We define a number of stochastic trend formulae (e.g. weakly increasing, strictly increasing, weakly decreasing, etc.) and use them to extend the set of state formulae of Continuous Stochastic Logic. We show how stochastic trends can be implemented in a guarded-command style specification language for transition systems. We illustrate the application of stochastic trends with numerous small examples and then we analyse the AKAP model in order to characterise and show causality and pulsating behaviours in this biochemical system

Process algebra modelling styles for biomolecular processes

We investigate how biomolecular processes are modelled in process algebras, focussing on chemical reactions. We consider various modelling styles and how design decisions made in the definition of the process algebra have an impact on how a modelling style can be applied. Our goal is to highlight the often implicit choices that modellers make in choosing a formalism, and illustrate, through the use of examples, how this can affect expressability as well as the type and complexity of the analysis that can be performed

Model Revision from Temporal Logic Properties in Computational Systems Biology

International audienceSystems biologists build models of bio-molecular processes from knowledge acquired both at the gene and protein levels, and at the phenotype level through experiments done in wildlife and mutated organisms. In this chapter, we present qualitative and quantitative logic learning tools, and illustrate how they can be useful to the modeler. We focus on biochemical reaction models written in the Systems Biology Markup Language SBML, and interpreted in the Biochemical Abstract Machine BIOCHAM. We first present a model revision algorithm for inferring reaction rules from biological properties expressed in temporal logic. Then we discuss the representations of kinetic models with ordinary differential equations (ODEs) and with stochastic logic programs (SLPs), and describe a parameter search algorithm for finding parameter values satisfying quantitative temporal properties. These methods are illustrated by a simple model of the cell cycle control, and by an application to the modelling of the conditions of synchronization in period of the cell cycle by the circadian cycle

Biochemical Reaction Rules with Constraints

International audienceWe propose React(C), an expressive programming language for stochastic modeling and simulation in systems biology, that is based on biochemical reactions with constraints. We prove that React(C) can express the stochastic pi-calculus, in contrast to previous rule-based programming languages, and further illustrate the high expressiveness of React(C). We present a stochastic simulator for React(C) independently of the choice of the constraint language C. Our simulator must decide for a given reaction rule whether it can be applied to the current biochemical solution. We show that this decision problem is NP-complete for arbitrary constraint systems C, and that it can be solved in polynomial time for rules of bounded arity. In practice, we propose to solve this problem by constraint programming

A High-Level Language for Rule-Based Modelling

Rule-based languages such as Kappa excel in their support for handling the combinatorial complexities prevalent in many biological systems, including signalling pathways. But Kappa provides little structure for organising rules, and large models can therefore be hard to read and maintain. This paper introduces a high-level, modular extension of Kappa called LBS-κ. We demonstrate the constructs of the language through examples and three case studies: a chemotaxis switch ring, a MAPK cascade, and an insulin signalling pathway. We then provide a formal definition of LBS-κ through an abstract syntax and a translation to plain Kappa. The translation is implemented in a compiler tool which is available as a web application. We finally demonstrate how to increase the expressivity of LBS-κ through embedded scripts in a general-purpose programming language, a technique which we view as generally applicable to other domain specific languages

On the spontaneous stochastic dynamics of a single gene: complexity of the molecular interplay at the promoter

International audienceBACKGROUND: Gene promoters can be in various epigenetic states and undergo interactions with many molecules in a highly transient, probabilistic and combinatorial way, resulting in a complex global dynamics as observed experimentally. However, models of stochastic gene expression commonly consider promoter activity as a two-state on/off system. We consider here a model of single-gene stochastic expression that can represent arbitrary prokaryotic or eukaryotic promoters, based on the combinatorial interplay between molecules and epigenetic factors, including energy-dependent remodeling and enzymatic activities. RESULTS: We show that, considering the mere molecular interplay at the promoter, a single-gene can demonstrate an elaborate spontaneous stochastic activity (eg. multi-periodic multi-relaxation dynamics), similar to what is known to occur at the gene-network level. Characterizing this generic model with indicators of dynamic and steady-state properties (including power spectra and distributions), we reveal the potential activity of any promoter and its influence on gene expression. In particular, we can reproduce, based on biologically relevant mechanisms, the strongly periodic patterns of promoter occupancy by transcription factors (TF) and chromatin remodeling as observed experimentally on eukaryotic promoters. Moreover, we link several of its characteristics to properties of the underlying biochemical system. The model can also be used to identify behaviors of interest (eg. stochasticity induced by high TF concentration) on minimal systems and to test their relevance in larger and more realistic systems. We finally show that TF concentrations can regulate many aspects of the stochastic activity with a considerable flexibility and complexity. CONCLUSIONS: This tight promoter-mediated control of stochasticity may constitute a powerful asset for the cell. Remarkably, a strongly periodic activity that demonstrates a complex TF concentration-dependent control is obtained when molecular interactions have typical characteristics observed on eukaryotic promoters (high mobility, functional redundancy, many alternate states/pathways). We also show that this regime results in a direct and indirect energetic cost. Finally, this model can constitute a framework for unifying various experimental approaches. Collectively, our results show that a gene - the basic building block of complex regulatory networks - can itself demonstrate a significantly complex behavior

A Novel Method to Verify Multilevel Computational Models of Biological Systems Using Multiscale Spatio-Temporal Meta Model Checking

Insights gained from multilevel computational models of biological systems can be translated into real-life applications only if the model correctness has been verified first. One of the most frequently employed in silico techniques for computational model verification is model checking. Traditional model checking approaches only consider the evolution of numeric values, such as concentrations, over time and are appropriate for computational models of small scale systems (e.g. intracellular networks). However for gaining a systems level understanding of how biological organisms function it is essential to consider more complex large scale biological systems (e.g. organs). Verifying computational models of such systems requires capturing both how numeric values and properties of (emergent) spatial structures (e.g. area of multicellular population) change over time and across multiple levels of organization, which are not considered by existing model checking approaches. To address this limitation we have developed a novel approximate probabilistic multiscale spatio-temporal meta model checking methodology for verifying multilevel computational models relative to specifications describing the desired/expected system behaviour. The methodology is generic and supports computational models encoded using various high-level modelling formalisms because it is defined relative to time series data and not the models used to generate it. In addition, the methodology can be automatically adapted to case study specific types of spatial structures and properties using the spatio-temporal meta model checking concept. To automate the computational model verification process we have implemented the model checking approach in the software tool Mule (http://mule.modelchecking.org). Its applicability is illustrated against four systems biology computational models previously published in the literature encoding the rat cardiovascular system dynamics, the uterine contractions of labour, the Xenopus laevis cell cycle and the acute inflammation of the gut and lung. Our methodology and software will enable computational biologists to efficiently develop reliable multilevel computational models of biological systems